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ABSTRACT 



Aims. Only for very few f5 Cephei stars has the behaviour of the magnetic field been studied over the rotation cycle. During the 
past two years we have obtained multi-epoch polarimetric spectra of the fj Cephei star V1449 Aql with SOFIN at the Nordic Optical 
Telescope to search for a rotation period and to constrain the geometry of the magnetic field. 

Methods. The mean longitudinal magnetic field is measured at 13 different epochs. The new measurements, together with the previous 
FORS 1 measurements, have been used for the frequency analysis and the characterization of the magnetic field. 
Results. V 1449 Aql so far possesses the strongest longitudinal magnetic field of up to 700 G among the p Cephei stars. The resulting 
periodogram displays three dominant peaks with the highest peak at / = 0.0720 d"' corresponding to a period P = 13'!893. The 
magnetic field geometry can likely be described by a centred dipole with a polar magnetic field strength B^ around 3 kG and an 
inclination angle fi of the magnetic axis to the rotation axis of 76±4°. As of today, the strongest longitudinal magnetic fields are 
detected in theySCephei stars V1449Aql and ^' CMa with large radial velocity amplitudes. Their peak-to-peak amplitudes reach 
~90kms"' and ~33kms"', respectively. Concluding, we briefly discuss the position of the currently known eight magnetic /? Cephei 
and candidate p Cephei stars in the Hertzsprung-Russell (H-R) diagram. 

Key words, stars: early-type — stars: fundamental parameters - stars: individual: V1449 Aql - stars: magnetic field - stars: oscilla- 
tions — stars: variables: general 



1. Introduction 

■Before the CoRoT mission, the B1.5 II-III star V 1449 Aql 
' (=HD 180642) was known as a /? Cephei pulsator with a high- 
amplitude radial mode of frequency 5.487 d"' (Waelkens et al. 
' I1998I I. White light photometry from space revealed the rich 
, frequency spectrum of this target (Degroote et al. I2009I I with 
the first detection of solar-like oscillations in a massive star 
(BeUcacem et al. 2009). Additional ground-based, multi-colour 
photometry and high-resolution spectroscopy were used to put 
' constraints on the atmospheric parameters and the chemical 
' abundances, and to identify three modes of low amplitudes 
(Briquet et al. I2009I I. Notably, the dominant mode does not 
behave sinusoidally, and the first velocity moment (y') has a 
peak-to-peak amplitude of about 90 km/s. This is the third 
largest peak-to-peak radial-velocity amplitude measured among 
yS Cephei stars, after BW Vul (e.g. Crowe & Gillet [T989] l and 
cr Sco (Mathias et al. 1991 ). A detailed NLTE abundance anal- 
ysis showed a nitrogen excess of ~0.3 dex with respect to the 
average nitrogen abundance in B-type stars in the solar neigh- 
bourhood. Previous abundance studies of hot stars indicate that 
the appearance of such an excess may be linked to the presence 



* Based on observations obtained at the 2.56-m Nordic Optical 
Telescope on La Palma and ESO Prgs. 077.D-031 1 and 178.D-0361. 



of a magneti c field (e.g.. Morel et al. l2006l l2008l Przybilla & 
Nieva BUTOt . 

The first spectropolarimetric observations of V 1449 Aql 
were obtained with the multi-mode instrument FORS 1 at the 
VLT on two consecutive nights in 2007 and revealed a change in 
polarity from one night to the next and the presence of a weak 
magnetic field on the second night of observations (Hubrig et al. 
2009). In the present work we discuss 13 new spectropolarimet- 
ric observations acquired during the last two years at the 2.56 m 
Nordic Optical Telescope using the SOFIN echelle spectrograph. 
The new longitudinal magnetic field measurements are used to 
determine the stellar rotation period and to put constraints on the 
magnetic field geometry. 



2. Magnetic field measurements and period 
determination 

Multi-epoch series of polarimetric spectra of V1449Aql with 
S/N >200 were obtained with the low-resolution camera (R = 
A/AA X 30000) of the high-resolution echelle spectrograph 
SOFIN (Tuominen et al. 119991 ) mounted at the Cassegrain fo- 
cus of the Nordic Optical Telescope (NOT) in 2009 August- 
September and in 2010 July. We used a 2K Loral CCD detector 
to register 40 echelle orders partially covering the range from 
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Table 1. Magnetic field measurements of VI 449 Aql with 
SOFIN. 



JD- 


Phase 


<s,> 


JD 


Phase 


<fiz> 


2450000 




[G] 


2450000 




[G] 


4343.659 


0.790 


-55±33 


5397.567 


0.648 


-430±89 


4344.584 


0.856 


166±41 


5398.530 


0.717 


-255±211 


5071.396 


0.171 


-21±95 


5400.574 


0.864 


168+168 


5073.391 


0.314 


-380±159 


5402.577 


0.008 


421±115 


5076.389 


0.530 


-647±153 


5405.646 


0.229 


-204±131 


5078.429 


0.677 


-326±171 


5407.559 


0.367 


-416±133 


5080.392 


0.818 


-115±102 


5409.648 


0.517 


-663±427 


5395.564 


0.504 


-729±232 









Notes. In the first two lines we list the earlier FORS 1 measurements 
published by Hubrig et al. C2009j . Phases are calculated according to 
the ephemeris of JD = 2455305.21 + 13.893 E. All quoted errors are Icr 
uncertainties. 




Fig. 1. The periodogram of V 1449 Aql built from the longi- 
tudinal magnetic field measurements with the most prominent 
frequency at / = 0.0720 d"', corresponding to a period of 
13.893 days. The lower thin line presents the window function. 



3500 to 10000 A with a length of the spectral orders of about 
140 A at 5500 A. Two exposures with the quarter- wave plate an- 
gles separated by 90° are necessary to derive circularly polarized 
spectra. The spectra are reduced with the 4A software package 
(Ilvinl2000). Wavelength shifts between the right- and left-hand 
side circularly polarized spectra were interpreted in terms of a 
longitudinal magnetic field {B^), using the moment technique 
described by Mathys ( 1 994 J . 

The magnetic field measurements with SOFIN and FORS 1 
are presented in Table [1] In the first column we list the heliocen- 
tric Julian dates for the middle of the spectropolarimetric obser- 
vations. The phases of the measurements of the magnetic field 
are shown in Col. 2, and in Col. 3 we present the mean longitu- 
dinal magnetic field {B^). An example of Stokes / and Stokes V 
profiles of N II in the wavelength region AA 4775^800 is pre- 
sented in Fig. lA.ll 

The frequency analysis was performed on the 13 SOFIN lon- 
gitudinal magnetic field measurements and two previous FORS 1 
measurements using a non-linear least-squares fit of the multiple 
harmonics utilising the Levenberg-Marquardt method (Press et 
al. 119921) with an optional possibility of pre-whitening the trial 
harmonics. To detect the most probable period, we calculated the 
frequency spectrum and for each trial frequency we performed 
a statistical F-test of the null hypothesis for the absence of peri- 
odicity (Seber [T977l l. The resulting periodogram displays three 
dominant peaks with the highest peak at / = 0.0720 d' corre- 
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Fig. 2. Phase diagram with the best sinusoidal fit for the longi- 
tudinal magnetic field measurements with a reduced ;\f^ - 0.46. 
The residuals (Observed - Calculated) are shown in the lower 
panel. The two FORS 1 measurements appear in green. 

sponding to a period P - 13'?893. In the framework of the rigid 
rotator model usually assumed for magnetic stars on the upper 
main-sequence, the period of the magnetic field variation corre- 
sponds to the stellar rotation period. In Fig. [T] we show that two 
lower peaks appear at / = 0.069162d-i and / = 0.059689d-i 
with equivalent periods of P = 14^459 and P - 16'?753, respec- 
tively. Since two measurements in our data set were obtained 
with the low-resolution FORS 1 spectrograph, we carried out 
a few tests to evaluate the significance of these measurements 
for the period determination. We find that omitting the FORS 1 
measurements from the period search leads to the dominance 
of a frequency peak coiTesponding to a period of 14'?496. On 
the other hand, the period search without weighting tied to the 
measurement accuracies leads to a much stronger dominance of 
the frequency peak / - 0.0720 d' with the equivalent period of 
P - 1 3'?893 . The referee of this paper informed us that his imple- 
mentation of a CLEAN analysis leads to a detection of the same 
periods, but favors the period of P = 14'?459. These results show 
that the periodicity determination may depend on the technique 
used and that the relevance of the frequencies / = 0.069 1 62 d' 
and / - 0.059689 d"' is difficult to ascertain at the present stage 
without additional magnetic data. The derived ephemeris for the 
detected period corresponding to the highest frequency peak is 

^g^^pos.extr. ^ HJD2455305.21 + 0.16 + 13.893 + 0.003£. (1) 

The phase diagram of the magnetic field measurements for the 
determined period is presented in Fig. |2] 

The measurements of the longitudinal magnetic field in 
V 1449 Aql require, however, special care in treating the data, 
as the timescale of a complete single observation of the length 
of ~40 min with SOFIN, including two subexposures, is compa- 
rable to the timescale of the pulsation variability. As we mention 
above, two subexposures with the quarter-wave plate angles sep- 
arated by 90° are needed to derive circularly polarized spectra. 
To remove the instrumental effects of small misalignments, dif- 
ferences in transmission, etc., the complete observation of a star 
usually consists of at least two successive subexposures with the 
waveplate rotated by 90° to exchange the positions of the two 
spectra on the CCD. The final Stokes / and V spectra are then 
calculated from the polarized spectra by taking the respective av- 
erages. To obtain Stokes /, the two average corresponding spec- 
tra, in our case the average of (/ + V)o and (/ -i- V')9o are added 
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Fig. 3. Variability of the output spectra in two SOFIN subexpo- 
sures taken with the quarter-wave plate angles separated by 90° 
taken ai'ound HJD 2455398.530. The lower two spectra, (/ + V)o 
and (/- y)o, correspond to the first subexposure, while the upper 
spectra, (/ - V)go and (/ + y)9o, correspond to the second subex- 
posure. The strong effect of pulsations on the line profile shapes 
and the line positions is clearly visible between the spectra of 
the first subexposure with a duration of 20 min and the spectra 
of the second subexposure with the same duration. 
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Fig. 5. Velocity moments of line profiles: equivalent width (mO), 
radial velocity (ml), line width (m2), asymmetry (m3), and kur- 
tosis (m4). In each panel five curves are plotted corresponding to 
five spectral lines of various elements: He 1 471 3. 1, C II 5133.1, 
N II 5667.6, O II 4592.0, and Si III 4553.6, shown from bottom 
to top. 



Table 2. Magnetic field model for V1449 Aql. 



<Sz> 


[G] 


-167.5±27.3 


^<-B.> 


[G] 


498.1±37.1 


P 


[d] 


13.893±0.003 


R 


[^o] 


11.2±1.4 


Veq 


[kms-'] 


41±5 


V Sin I 


[kms-'] 


24±2 


i 


n 


36±6 


P 


n 


76+4 


Bi 


[G] 


3100±300 



Fig. 4. Time series of FEROS spectra showing pulsational line 
profile variability in the spectral region 4590-4615 A. The pul- 
sation phase zero is at the bottom. 

with the average of (/ - V)q and (/ - y)9o, while Stokes V is 
calculated from the difference of the two spectra, the average of 
{I+V){) and (/-i-y)9o and of (/-y)o and (/-y)9(). However, owing 
to very strong changes in the line profile positions and shapes in 
the strongly pulsating VI 449 Aql, the method of using respec- 
tive average spectra would lead to erroneous wavelength shifts 
and thus to wrong values of the longitudinal magnetic field. The 
behaviour of all four spectra in both subexposures taken around 
HJD 2455398.530 is illustrated in Fig.|3] To remedy this situa- 
tion, we decided to carry out the measurements of the line shifts 
between the spectra (/ + V)q and (/ - V)q in the first subexposure 
and the spectra (/ - y)9o and (/ + V)c)q in the second subexposure 
separately. The final longitudinal magnetic field value was then 
calculated as an average of the measurements for each subexpo- 
sure. 

The line profiles undergo very significant pulsational vari- 
ability, and any rotational modulation of the observed profiles is 
completely masked by the pulsational modulation. In Fig. IA.2I 
we demonstrate the behaviour of the O ii and He i line profiles 
in SOFIN spectra in the spectral region 4695^7 15 A over the 
rotation cycle. The inspection of 224 FEROS spectra used in the 



study of Briquet et al. (I2009I I reveals that all lines behave in a 
rather similar way. In Fig.|4]we present a clear pulsation pattern 
in the spectral region 4590^615 A. The moment variations for 
five lines belonging to different elements, namely He I 4713.1, 
CII 5133.1, Nil 5667.6, Oil 4592.0, and Silll 4553.6, shown 
in Fig. |5] however, exhibit slight diff'erences between the ele- 
ments. The equivalent widths of O II and Si III are larger at the 
phase of minimum radius, while C II and N II are larger at the 
maximum radius. The equivalent width of He I is almost con- 
stant. At minimum radius the elements N, C, and O show line 
profiles narrower than those of Si and He. Also the amplitude of 
the radial velocity curve of the He I line is about 4 km s"' greater 
than for the C II line and 2 km s ' more than for the N II line. 

The longitudinal field of VI 449 Aql can be adequately rep- 
resented by a single-wave variation curve during the stellar ro- 
tation cycle (see Fig. |2]i, indicating a dominant dipolar contri- 
bution to the magnetic field topology. Assuming that V1449 Aql 
is an oblique dipole rotator, the magnetic dipole axis tilt p is 
consti-ained by (BJ™" / (Bz)™" ^ cos(/ + y3)/ cos(/ - jS). In 
Table |2] we show the mean value {Bj} and the semi-amplitude 
of the field variation A(By) for V1449Aql in rows 1 and 2. 
Values for Vt - 44km s"', which is the total line broadening, 
and Veq = 38 + 15 kms"' were published by Briquet et al. 
(2009'). The stellar radius, R = 11.2 ± \ARq, was calculated 
using main-sequence evolutionary CLES models (Scuflaire et 
al. [2008 ). Obviously, the accuracy of the determination of Veq 
from the best-fit-solutions of the spectroscopic mode identifica- 
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tion described in the work of Briquet et al. (I2009I I is rather low. 
However, the equatorial rotation velocity can also be calculated 
using the relation Veq = 50.6 RIP, where R - \\ .2 + \ A Rq \s 
the stellar radius in solar units and P - 13.893 d the period in 
days. Using this formula, we obtain Vgq = 41 + 5 km s"' . On the 
other hand, Lefever (private communication) obtained v sin / = 
24kms"\ using the different procedure described in Lefever et 
al. (2010), which takes macroturbulent broadening into account. 
Our own estimate of the vsin; value, vsin; = 23.5 km s"' us- 
ing the technique of Diaz et al. (120101 ) fully agrees with that of 
Lefever The combination of vsin; - 24kms~' with a typical 
estimation error of 2 km s"' and v. 



eq 



41+5 km s" leads to the 
inclination angle / = 36 + 6°. In the last two rows, we list the 
parameters of the magnetic field dipole model. 

3. Discussion 

Using FORS 1/2 and SOFIN longitudinal magnetic field mea- 
surements collected in our recent studies, we were able to detect 
a weak mean longitudinal magnetic field of a few hundred Gauss 
in fourySCephei stars, JCet, ^' CMa, 15CMa, and V1449Aql, 
and in two candidate j8 Cephei stars, a Pyx and e Lup (Hubrig et 
al. 120061 120091 |20TTI) . Before we started our systematic search 
for magnetic fields in pulsating B-type stars, a weak magnetic 
field was detected in two other ySCephei stars, in the proto- 
type of the class, /3Cep itself, by Henrichs et al. (.2000) and 
in V2052Oph by Neiner et al. (120031 ). The strongest longitudi- 
nal magnetic fields were detected in the two /? Cephei stars with 
the large radial velocity amplitudes, namely in V 1449 Aql and 
^' CMa, for which peak-to-peak amplitudes reach ~90kms"' 
and ~33kms ', respectively. The rather strong magnetic field 
in the atmosphere of these two stars and the rather low v sin / 
values, certainly make these stars the most suitable targets for 
studying the impact of magnetic fields on stellar rotation, pulsa- 
tions, element diffusion, and frequency patterns. 

In the frequency analysis of the CoRoT data obtained for 
V1449 Aql (Degroote et al. l2009b . the frequencv f = 0.068790+ 
0.000393 d"' is, within the errors, close to the frequency / = 
0.0691 62 d"' detected in our study of the periodicity using mag- 
netic field measurements. Additional magnetic field measure- 
ments will be worthwhile for proving whether this photomet- 
ric frequency corresponds to the rotation/magnetic period of this 
star. 

Our studies indicate that dipole models provide a satisfac- 
tory fit to the magnetic data and among the presently known 
magnetic jSCephei stars, f ' CMa and V1449Aql, possess 
the largest magnetic fields, with a dipole strength of several 
kG. The position of the currently known eight magnetic 
jS Cephei and candidate /? Cephei stars in the logTgff-logg 
diagram is presented in Fig. |6] together with the boundaries 
of the theoretical instability strips calculated for differ- 
ent metallicities (Z - 0.01 and Z = 0.02) and using the 
OP opacities (http://cdsweb.u-strasbg.fr/topbase/op.html, 
see also Miglio et al. .2007) and OPAL opacities 
(http://opaiopacity.linl.gov/opal.htmI). The stellar funda- 
mental parameters and their literature sources are compiled in 
Table lA.ll We give in the same table the information on the 
dominant mode and the radial velocity amplitude, if known. 
The extent and the position of the instability strip are highly 
dependent on the opacities used and on the metallicity adopted. 
Recent studies indicate, however, a lower metallicity than 
Z = 0.02 for early B-type stars in the solar neighbourhood with 
Z = 0.014 (e.g. Przybilla et al. l2008i) . Although the sample of 
magnetic pulsating stars is still very small, and it is difficult to 
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Fig. 6. The position of the fi Cephei and candidate y6 Cephei stars 
in the H-R diagram. Upper panel: The boundaries of the the- 
oretical instability strips are calculated using the OP opacities. 
Lower panel: Boundaries are calculated using OPAL opacities. 
Full lines correspond to strips for metallicity Z = 0.01 and dotted 
lines to strips with metallicity Z = 0.02. The thick and thin lines 
correspond to the boundaries of the (3 Cephei and SPB instability 
regions, respectively. 



determine any trend in the distribution of the stars in the H-R 
diagram, it might seem that magnetic stars are clustered rather 
close to the instability strip boundaries (full lines in left panel of 
Fig. |6]). Certainly, additional future magnetic field studies of a 
larger sample of yS Cephei stars are needed to study the indicated 
trend. 
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Fig.A.2. Line profile variabiHty of On ^^4699, 4705, 4710, 
and He I A4-113 over the rotation cycle on SOFIN spectra. The 
spectra are presented with phases increasing from bottom to top 
and offset in the vertical direction for clarity. 



Appendix A: Additional material 
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Table A.l. Fundamental parameters of jSCephei and candidate ySCephei stars with detected magnetic fields. 

Object Teff log g Reference dominant mode RV amplitude Reference 

name [K] [kms"'] 



(5Cet 


23000±1000 


3.90+0.10 


1 


/■= 6.2026 d-' 


7.44 


5 


^' CMa 


27000±1000 


3.80+0.10 


1 


/ = 4.77153 d-' 


16.4 


6 


15CMa 


24000+1000 


3.40+0.10 


1 


/ = 5.42d-' 




7,8 


ffPyx 


22900+300 


3.60+0.05 


2 


/3 Cep candidate 




9 


eLup 


19300+900 


3.89+0.20 


3 


P Cep candidate 




9 


V2052Oph 


23000+1000 


4.00+0.20 


1 


/ = 7.148 d-' 


6.7 


10 


V 1449 Aql 


24500+1000 


3.45+0.15 


4 


/ = 5.48694d-' 


44.0 


4 


/3Cep 


25000+1000 


3.80+0.10 


1 


/ = 5.250 d-' 


21 


11 



No tes. 1 - Lefever et al. l l20T0t, 2 - Przybilla et al. ( [20T0b , 3 -H ubrig et al. ( [2009] >, 4 -Br iquet et al. ( |2009) , 5 - Ae rts et al. ( fT992t . 6 - Sa esen e t 
al. f2006] l, 7 - Heynderickx fl9^ , 8 - Heynderickx et al. ( fT994l . 9 - Telling et al. ( l2006l . 10 - Neiner et al. ( 12003) , and 1 1 - Aerts et al. ( fT994t 



